Spatial Dependence of Entropy in Trapped Ultracold Bose Gases 
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We find a new physical regime in the trapped Bose-Hubbard Hamiltonian using time-evolving 
block decimation. Between Mott-insulating and superfluid phases, the latter induced by trap com- 
pression, a spatially self-organized state appears in which non-local entropy signals entanglement 
between spatially distant superfluid shells. We suggest a linear rather than harmonic potential as 
an ideal way to observe such a self-organized system. We also explore both quantum information 
and thermal entropies in the superfluid regime, finding that while the former follows the density 
closely the latter can be strongly manipulated with the mean field. 
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Entropy is a fundamental concept in physics, chem- 
istry, and information theory. It is as foundational as en- 
ergy for many-body systems, whether classical or quan- 
tum. For example, consideration of entropy leads to the 
second law of thermodynamics, which is necessary to un- 
derstand the efficiency of engines; thus entropy has had 
practical applications from the time of its invention in the 
nineteenth century. In ultracold quantum gases, an ex- 
cellent playground for quantum many-body physics and 
therefore basic many-body concepts such as entropy, the 
relationship between entropy and temperature for Bose 
and Fermi statistics is central to the concept of adia- 
batic cooling [1]; the latter led to the discovery of the 
atomic Fermi superfluid and profound new insight into 
how fermions pair to become bosons [l, 0, 0]- Moreover, 
entropy serves as an entanglement measure in such quan- 
tum systems [5]. However, in all these cases entropy is 
treated as a scalar - a bulk measure. In this Letter, we 
show that taking entropy as spatially dependent, whether 
as a zero temperature quantum information (QI) con- 
cept or a finite-temperature statistical mechanical con- 
cept, leads to new applications and new understanding 
of Bose gases. 

Specifically, we start with the trapped Bose-Hubbard 
Hamiltonian (TBHH) [6] at zero temperature. Us- 
ing time-evolving block decimation (TEBD), we explore 
three varieties of QI entropy: local von Neumann entropy, 
also called entropy of entanglement; block entropy; and 
mutual entropy, which is non-local. The TBHH has two 
axes in the thermodynamic limit, namely chemical poten- 
tial and tunneling, both in ratio to interaction. Consid- 
eration of the local entropy shows that a third axis, trap- 
ping energy in ratio to interaction, is necessary to under- 
stand the TBHH. We find a spatially self-organized state 
in this fundamentally finite system. The characteristic 
feature of this new physical regime intermediate between 
superfluid (SF) and Mott Insulator (MI) is high non-local 
entropy between SF shells, as shown in Fig. [D^c). Al- 
though this effect appears already for the harmonic po- 
tential, we suggest the linear potential as an ideal place in 
which to explore this idea. Then we go on to explore the 
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FIG. 1: (color online) Mutual entropy in the "wedding cake^' 
regime of harmonic trap, where Mott-insulating layers are sep- 
arated by superfluid shells, (a) The local entropy (dashed red 
curve) is oscillatory but non-zero in the MI layers; the density 
is shown for comparison (solid blue curve), (b) The mutual 
entropy is oscillatory along the diagonal, indicating mainly 
nearest- neighbor entropy; the two bright spots are SF regions, 
(c) When the system is compressed by the trap, highly non- 
local entropy appears, (d) Further compression can lead to 
nonlocal entropy being dominant and completely SF system. 



SF regime. For a Thomas-Fermi (TF) profile in the 
lattice plus harmonic trap, all three QI entropies closely 
follow the density. On the other hand, thermal entropy 
behaves completely differently. In the local density ap- 
proximation (LDA) the entropy falls away rapidly in the 
region outside the condensate, but penetrates deeply into 
the condensate, having a small attrition of typically only 
about 50% in the center; thus a Feshbach resonance can 
be used to strengthen the mean field and push the en- 
tropy out to the edges of the condensate, thereby max- 
imizing the efficiency of evaporative cooling below the 
critical temperature. 
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The first part of our study explores the TBHH outside 
the LDA. Quantum phase transitions in finite systems 
where the LDA fails is key not only to fundamental fea- 
tures of quantum gases such as superfluidity Q but also 
in nuclear physics [si]. For the latter part of our study 
concerning thermal entropy in the LDA, a very recent 
work by Bernier et al. [10] also proposes spatial manip- 
ulation of entropy, although purely in the LDA and for 
fermions. They predict that an order of magnitude tem- 
perature decrease can be achieved, their main tool being 
dynamical mean field theory in the LDA. For bosons, the 
study of thermal entropy as a bulk quantity has a long 
history and is well summarized in Dalfovo's review 0. 

We begin with a brief discussion of QI entropies. The 
direct quantum analog of the thermal, or Gibbs entropy 
is the von Neumann entropy, also called the entropy of 
entanglement: 5'vn = — Trplog/3, where p is the state 
matrix, and indeed, for finite temperature we know that 
'S'cibbs in the thermodynamic limit. However, at 
zero temperature the von Neumann entropy is zero since 
the system is in a pure state. Instead, we can ask about 
the local entropy based on local measurements of the on- 
site density matrix pi = Tij-Lip in a lattice model. In 
this vein, we define three distinct QI entropies: 

^local ^ _Trp,l0gPi, (1) 

^biock ^ _Trp<,log/5<,, (2) 

cmutual Q olocal olocal /Q^ 

^ij — ^ij ^3 •> 

where Sij means tracing over all sites but two, and p<i 
means tracing over all sites to the right of site i. The 
local entropy is the nearest QI equivalent to the thermal 
spatial entropy density S^^{x). The block entropy results 
from a bipartite splitting of the lattice made at site i. 
The mutual entropy or mutual information [5] assumes a 
measurement made at two sites only. 

We now proceed to consider a harmonically trapped 
Bose gas in a lattice. The TBHII, which assumes the 
tight-binding lowest-band approximation, is given by 

-^VoE^^''rl^,n^, (4) 

where J is the hopping strength, U is the interaction 
strength, Voi^ is a power law potential with p G |Z|, and 
hi, b] are bosonic destruction/creation operators satisfy- 
ing bosonic commutation relations on-site and commut- 
ing between sites. We use TEBD in ID to study the 
three QI entropies of Eqs. ([I])-®. TEBD has two con- 
vergence parameters, the on-site dimension and the 
number of eigenvalues x retained in the reduced density 
matrix. In all cases we have checked convergence by at 
least doubling both d and x and seeing if we obtain the 
same results. Our values of d range from 8 to 12, meaning 
7 to 11 bosons per site (8 to 12 including the vacuum), our 
values of x range from 30 to 120, our number of sites is 
from L = 30 to 100, we use a number-conserving version 
of TEBD to make these large (d,x)-values tractable, and 



we implement imaginary time propagation in the Suzuki- 
Trotter scheme to obtain the ground state. 

We first consider the MI regime, in which it has been 
shown that the harmonic trap causes the system to spa- 
tially traverse the TBHH phase diagram in the LDA, 
leading to a layered structure with SF shells separating 
MI layers [l2|, ; this has been dubbed "the wedding 
cake." Thermal entropy has been shown to be concen- 
trated in the SF shells in this case, since the MI layers 
are a near-product state, and therefore have very small 
entropy; such considerations have led to new cooling 
schemes [l^. In contrast, for QI entropy, if the MI layer 
is less than about 15 sites, the entropy is much larger even 
far into the Mott lobe (more correctly, a "Mott claw" in 
ID), due to finite size effects [l5|; in this region, which is 
quite typical in experiments, the assumption of the LDA 
is of limited value. However, we still find that the lo- 
cal entropy is large in SF shells and small in MI layers, 
varying by a factor of between 2 and 3 for J/U = 1/16, 
Vq/U = 0.0125, and = L = 30 bosons, for example, 
as shown in Fig.[T](a). This is qualitatively similar to the 
thermal entropy result. What about nonlocal entropy? 
In Fig. [Hb) we show the mutual entropy corresponding 
to the parameters of Fig. (TJa). The wedding cake struc- 
ture is apparent: the mutual entropy is mainly nearest- 
neighbor and small in the MI layers, while it is larger and 
slightly more delocalized in the SF shells. 

We have so far taken the harmonic trap to be suffi- 
ciently weak so that, even if not in the LDA, the MI layers 
are still many sites long. What happens when we com- 
press the trap? This compression requires a third axis 
in the TBHH phase diagram and is a natural direction 
to explore in experiments. For sufficiently strong traps 
the incompressibility of the MI phase is overcome and 
the system melts into a SF. Considering such compres- 
sion as a continuous process we observe in simulations 
that the MI layers melt gradually, shortening from many 
sites to just a few. When each MI layer is reduced 2 or 3 
sites we find a completely different form for the mutual 
entropy, as shown in Fig. [Tfc) {Vq/U = 0.025) and (d) 
{Vq/U = 0.5). The entropy becomes highly nonlocal due 
to penetration of SF tails through intervening MI layers. 

In the harmonic trap this phenomenon is common but 
not regular, because a small amount of SF shell is re- 
quired between each MI layer; due to the potential dif- 
ference between sites being non-uniform (harmonic) this 
does not always occur. In contrast, in a linear potential 
the alternation between SF shells and MI layers is more 
regular and generic, since special parameters are not re- 
quired and just compressing the trap naturally leads to 
this spatially self-organized state. We use the word "self- 
organization" in the sense that the system optimizes ar- 
rangement of narrow SF and MI layers. A typical ex- 
ample is shown in Fig. [2](a)-(b), where seven SF regions 
produce a nested square-like pattern in the mutual en- 
tropy. In general, when the MI layers are a few sites long 
then there is a strong non-local mutual entropy, indicat- 
ing that the MI layers act as effective potential barriers 
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FIG. 2: (color online) New physical regime of the Bose- 
Huhhard Hamiltonian for a linear potential, (a) Density and 
local entropy, and (b) mutual entropy - compare to the har- 
monic potential explored in Fig. [1] (c) Compressibility k, for 
J/U = 1/8 (red points, solid curve) and 1/16 (black points, 
dashed curve) as a function of the trap strength. The inner 
panel shows a zoom around the value of Vb/t/ = 0.24 used for 
panels (a) and (b). 



for the SF regions. Thus there are three physical regimes: 
(1) for a weak linear trap the system is a MI; (2) for a 
stronger potential the bosons form alternating SF shells 
and MI layers, with an average wavelength A which de- 
creases monotonically as Vq increases; and (3) completely 
melted, or SF, for very strong trap compression. In 
Fig. [21(c) we show the compressibility = —(j~^d(j/d{Vo) 
where a is obtained by a Gaussian fit. After the abrupt 
change between from pure MI to wedding cake structure, 
the onset of strong oscillations indicates regime (2); such 
oscillations become longer wavelength inVo/U and damp 
out in the purely SF regime (3). The sudden jumps in 
have been discussed previously in the context of changing 
interaction [l^; here we observe the same kind of effect 
caused by compressing with the trap. The jumps are due 
to new superfluid regions entering the system. Smaller 
J/U leads to a stiffer response and therefore more sudden 
changes in k,. The zoom shows that these changes are in 
fact smooth, simply more abrupt in some regions. The 
outer panel of Fig.[2](c) has a resolution of 500 points and 
the derivative for is approximated to second order. 

In order to make the comparison between QI entropy 
and the more familiar thermal entropy, we turn to the 
harmonic trap and the completely SF regime of Eq. (|4]), 
J/ZU > 1, where Z is the coordination number. We 
consider the TF regime in which the condensate takes a 
characteristic harmonic shape. In this weakly discretized 
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FIG. 3: (color online) Quantum information entropies in the 
superfluid regime in a ID lattice plus harmonic trap. Shown 
are (a) density, (b) local entropy, (c) block entropy, and (d) 
mutual entropy, all on 100 sites. Panels (a)-(c) show N — 
100, 200, 300 atoms (red dash-dotted, blue dashed, and solid 
black curves); panel (d) shows 300 bosons only. 



limit the lattice can provide a convenient discretization 
to approach the continuum problem in the finite trapped 
system by holding U /vJ = ^ fixed, where v is the filling 
and 7 is the Tonks parameter, and letting the number of 
sites be large [l^. We approach the SF continuum but 
remain still in the weakly discretized regime, with 100 
sites and hundreds of bosons, where the LDA is mean- 
ingful. Figure [3] shows the density and three QI entropies 
in this case, with J = /7 = 1 and = 0.004. The QI 
entropies closely follow the density - note the TF pro- 
file, indicating that we are in the near-continuum limit. 
Specifically, we find that the ratio of the full-width half- 
max (FWHM) of the entropies, local or non-local, to the 
FWHM of the density, is approximately constant over a 
wide range of parameters, studying J/U from 1/2 to 3 
and a variety of system sizes and particle numbers - the 
ratios of the FWHMs varies by only a few percent in any 
given case. 

As we show in the following, the thermal entropy be- 
haves completely differently from the QI entropy; it is 
even expelled from the center of the trap. This is in 
contradistinction to the MI regime, where both QI and 
thermal entropy showed the same behavior. The spatial 
entropy or local entropy density S^^{x) of a weakly in- 
teracting harmonically trapped Bose gas in the LDA in 
the semiclassical approximation and using the TF pro- 
file [l|, 0], i.e., assuming T <C Tc and 7 <C 1 for our 
quasi- ID system [l8|, is given by 

=kBj^ [f^e^i^ip, x) - In (1 - e"'^^-)] , (5) 

where i^{p^x) = [exp(;^eb) — 1]""^ is the Bose weighting 
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FIG. 4: (color online) Thermal entropy density for a BEC in 
a harmonic trap in the Thomas-Fermi regime. Outset: Spa- 
tial entropy density in the inner and outer regions of a Bose- 
Einstein condensate, as functions of 2pgn{x) '2P[iJ. — V{x)] 
(solid black curve) and f3[V{x) — /j] (dashed blue curve). Note 
the rapid fall-off in the outer trap region compared to the slow 
fall-off in the inner, condensate region. Inset: Distribution of 
the entropy (solid curves) for parameter range of Oberthaler 
experiments (green solid, dashed red curves). 



factor and the energy of Boguliubov excitations is 



/2m - 



>R. 



(6) 

In Eq. dl]), /i is the chemical potential, R = (2/i/mcc;^)^/^ 
is the TF radius, n{x) is the mean-field density, m is 
the atomic mass, giB = hu;±as is the contact interaction 
strength with the s-wave scattering length, and uo <C 
cjx, with u;± the transverse trapping frequency. The TF 
regime is reached for (i ^ hw^ and involves neglecting 
the quantum pressure in the mean field theory, giving 
the condensate an inverted parabolic shape. 

In Fig. m we show the numerical solution of Eq. ([5]) both 
inside and outside the TF radius, each integral being per- 
formed according to Eq. (|6|) , with the units of entropy per 
volume encapsulated m A = kB{V2/Ti'^){m/ . Al- 
though the peak value of the entropy is always reached 
on the ^-axis of Fig.[4l which occurs at the TF radius, the 
fall-off is not the same outside and inside the condensate. 
In the outer region S^^{x)/A falls away very quickly as a 
function of (3\y(x) — /i]. In contrast, in the inner region 
it falls away much more slowly as a function of 2Pgn{x). 
The extent of the fall-off is given by the peak condensate 
density as 2/3gn{0) = in the TF regime. For the 
experimental parameters of the Oberthaler group's ^^Rb 
experiments, (3fi ranges from 2 to 10, uj± = 27r x 400 Hz, 



= 27r X 20 Hz, = 5.29 x 10"^ cm, ksT c:^ 10 to 20 nK, 
and there are 1000 to 2000 bosons. Then 2(3gn{0) = 4 to 
20 and the entropy falls off by between 70% and 93% from 
the edge of the condensate to the ce nter, a s sketched in 
the inset of Fig.[4]in "trap units" = ^/h/mO. A Feshbach 
resonance can be used to increase this range greatly, lo- 
calizing the trap entropy even more strongly at the trap 
edges for more effective evaporative cooling below the 
critical temperature. After an evaporative cooling step 
the interactions can be decreased to avoid strong losses 
incurred by the Feshbach resonance. Adiabatic and cycli- 
cal transport of the entropy in this manner, in conjunc- 
tion with a careful study of loss effects pj|], provides a 
subject for future theoretical and experimental investiga- 
tions. 

In past studies of the wedding cake structure, phase 
coherence between SF shells was neither expected nor 
observed. In the new highly-compressed regime we sug- 
gest, an interference experiment performed by turning 
off the trap and letting the system expand should result 
in an observable periodic interference pattern on top of 
the smeared-out pattern corresponding to the MI. In re- 
peated experiments the position of the peaks should be 
stable due to phase coherence. Recent 3D LDA studies of 
Josephson oscillations between SF layers separated by MI 
plateaus suggest that such phenomena are not restricted 
to ID [20,[2l|. 

In conclusion, we explored spatial control of entropy 
in trapped Bose gases. We used quantum information 
entropy to find a new axis of the trapped Bose-Hubbard 
Hamiltonian outside of the LDA, namely the harmonic 
trap strength. Between MI and SF phases, the latter in- 
duced by trap compression, we found a new spatially self- 
organized state signaled by large non-local entropy; we 
suggested the linear potential as an ideal trap in which 
to observe this effect. Turning to the weak-lattice su- 
perfluid regime where LDA is valid, we found that QI 
and thermal entropies do not display the same pattern. 
Although QI entropy closely follows the boson density, 
thermal entropy can be localized at the trap edges, de- 
pending on the mean field strength; this leads to new 
cooling applications for the future. 
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